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Regenerators have been proposed to increase thermal cycle 
efficiency of gas turbine power plants which are being considered 
for installation in certain classes of United States Navy ships. 

These regenerators are to preheat intake air and are currently designed 
as compact heat exchangers. The recommended heat exchanger configur- 
ations are either counterflow or crossflow. This analysis considers 
only the counter flow configuration. 

Some of the structural problems associated with rectangular, com- 
pact heat exchanger design are discussed and certain structural design 
parameters are determined analytically. A digital computer program is 
presented for the analysis of the pressure effects in an idealized, 
two dimensional, rectangular heat exchanger. This program is written 
for a matrix Analysis of an idealized heat exchanger with a cross 
section analogous to a multi-story building frame. 

This investigation was conducted by Robert L. Corbett and David 
W. Stubbs during the period from December 1965 to May 1966, at the 
United States Naval Postgraduate School, Monterey, California. 
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TABLE OF SYMBOLS 



a Height of the heat exchanger module, in 
b Width of the heat exchanger module, in 
c Height of one flow channel in the heat exchanger, in 
D Plate flexural rigidity (EI/l-v^), lb-in^ 

r\ 

E Modulus of elasticity, lb/ in 

f Distance between fins, in 

H Axial force, lb 

I Moment of inertia, in^ 

El for the vertical side of the elemental frame, lb-in^ 

J2 El for the horizontal side of the elemental frame, lb-in 
K Stiffness influence coefficient 
L Length of the heat exchanger module, in 
M Moment, in- lb 

m Frame member end moments, in- lb 
p Frame member end forces, lb 

o 

Q Heat exchanger internal pressure, lb /in 

t i Thickness of the vertical side of the elemental frame, in 

t 2 Thickness of the horizontal side of the elemental frame, in 
t a Thickness of the heat exchanger closure bars, in 

tjj Thickness of the heat exchanger top plate, in 

tj Thickness of the heat exchanger fins, in 

tp Thickness of the heat exchanger interior plates, in 

V Shear force, lb 

Poisson* s ratio 

9 Frame joint rotations, radians 

<f Frame joint deflections, in 
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1 . 



INTRODUCTION 



Several current projects in industry and one at the United States 
Naval Postgraduate School are actively concerned with the performance 
of gas turbine regenerators* The impetus for these investigations stems 
from a proposal to equip certain types of United States Navy ships with 
a power plant which will reduce, significantly, the delay time for 
transition from cruising speed to flank speed. 

This desired flexibility in speed is planned to be achieved either 
by augmenting or replacing the current steam turbine installations with 
diesel or gas turbine power plants. In general, the efficiency of the 
diesel is superior to that of the pure gas turbine installations; but 
the advantage in weight and space for gas turbine engines is significant. 
These advantages, coupled with increased efficiency gained through the 
use of regenerators, make the gas turbine appealing for shipboard 
installations. Therefore the design, operation, and maintenance of 
gas turbine regenerators are important in the proposed plan for gas 
turbine shipboard installation. 

Regenerator or heat exchanger requirements for different types of 
ships may differ in size, shape, and weight. A basic heat exchanger 
may be considered to consist of various combinations of small rectangular 
modules stacked or joined together in a manner that allows the complete 
heat exhcanger to meet the heat transfer, geometry, and weight require- 
ments for a particular installation. The rectangular shape is a log- 
ical choice for stacking and combining of modules and is simple to 
fabricate. 

Present design requirements specify maximum pressures up to ten 
atmospheres, temperatures ranging from 500 to 1000 °F, and a service 
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life ranging from 20,000 to 60,000 hours. This latter requirement is 
severe, since structural integrity problems have occurred in similar 
heat exchangers after operating for much shorter periods of time. 

The Bureau of Ships has suggested that the following areas be 
investigated: 

1. Gas side fouling and air side salt precipitation 

2. Thermal stress behavior 

3. Thermal fatigue criteria 

4. Materials evaluation 

5. Header design and flow distribution 

6. Fabrication and joining techniques 

7. Advanced heat exchanger designs 

8. Engine evaluation of advanced heat exchangers 

This list indicates some of the areas which require investigation in 
order to define the requisite design criteria and to determine the order 
of reliability to be specified for a compact heat exchanger. 

Preliminary stress investigations of typical rectangular heat 
exchangers have been accomplished to date. Few complete analyses of 
the complex state of stress within the heat exchanger core or supporting 
frame work have been made and those reported are, in general, sufficient 
only for the justification of the preliminary design. These result 
in heavier and more costly heat exchanger design, and possibly in 
reduced reliability. The authors believe that, through a careful 
analytical study to determine certain structural design parameters, 
heat exchanger design may be determined more readily. 

The stresses in different parts of a heat exchanger are due 
primarily to the following: CO 
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1. The difference in static and dynamic fluid pressures and 
atmospheric pressure. 

2. The thermal expansion or contraction of component parts 
under: 

a. Steady state conditions 

b. Transient or thermal shock conditions 

3. Method of fabrication. 

The values of these stresses are difficult to determine not only 
because of the complex structure but also since the value of pressure 
and temperature is a function of location and time. 

The purpose of this investigation is to develop a basic method 
to study the effect of pressure on the design. Some of the problems 
relating to structural integrity are discussed and an attempt is made 
to isolate certain design parameters. The effect of some such para- 
meters is qualitatively studied through the use of the high speed 
digital computer. 

The investigation reported in this paper was conducted by the 
authors during the period from December 1965 to May 1966 at the United 
States Naval Postgraduate School, Monterey, California. Sincere 
appreciation is extended to the faculty and staff of the Aeronautical 
Engineering Department for their assistance. The authors are deeply 
indebted to Professor C. M. Smith for his invaluable encouragement and 
assistance throughout the project. The authors also wish to express 
their thanks to Professor C. H. Kahr for introducing them to and in- 
structing them in structural analysis by matrix methods. 
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2 . 



GENERAL 



Heat exchangers of the cylindrical or tube type design have been 
analyzed and have proven to be structurally satisfactory. On the other 
hand, compact-rectangular heat exchangers have not been analyzed exten- 
sively for their structural action. 

A number of analytical approaches to the structural problem were 
explored in this investigation of pressure effects. Some of the methods 
proved to be not entirely satisfactory and seveial of these are outlined 
in Appendix II. 

The method of analysis adopted was limited to the core and was 
independent of the flange. The core is important from a heat transfer 
point of view, and is undoubtedly the most costly. The flange is not 
included in this investigation but its effect on structural design is 
sufficiently important to warrant further consideration. 

In the analysis of the core, cross sections were considered to be 
similar to multi-story building frames with the closure bars forming 
main structural columns, the fins acting as partition walls or columns, 
and the interior plates acting as floors. Exploratory design parameters 
were determined by considering a rectangular frame without interior 
tie-bars or constraints. The effect of these parameters on the lattice- 
like frame composing the core cross section was then investigated. 
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3 . 



ANALYSIS 



For the purpose of this investigation, a rectangular core module 
for a counter flow heat exchanger was postulated. In this module the 
geometry and material properties were assumed to be constant over the 
length; the pressure was assumed to be constant in all channels. The 
cross section was divided into a pattern of small rectangular passages 
by the interior plates and fins. All end discontinuities resulting from 
attachment of the core to flanges were neglected. This idealized core 
module is presented in Figure 1. 
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FIGURE 1 

IDEALIZED COUNTERFLOW CORE MODULE 
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This module represents only a small part of the complete regenerator 
heat exchanger. Different configurations would have a number of these 
modules arranged in a manner suitable to the operational environment 
of the ship. 

Close inspection of the idealized core module reveals that the 
fundamental element involved in the cross section is an open rectangular 
frame. Many of these basic frames are contained in a module cross section 
and two or more modules are stacked together to form a complete heat 
exchanger . 

Since the fundamental element is an open rectangular frame, an 
analysis was performed on a simple frame in order to isolate parameters 
which may affect or define structural considerations for design. DU 
In this analysis the frame was assumed to be loaded with a uniform 
internal pressure and to have two axes of symmetry as shown in Figure 2. 
While these assumptions are not strictly applicable to the arrangement 
of elemental frames in the core module, they do facilitate the analysis 
and provide information as to which parameters affect the stresses in 
each member of the frame. 
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FIGURE 2 

ELEMENTAL FRAME, TYPICAL 
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The two axes of symmetry which are shown in Figure 2 permit the 
analysis to be made for one quarter of the elemental frame; the required 
boundary conditions are defined through considerations of symmetry. 

The frame representing one quarter of the elemental frame is shown in 
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FIGURE 3 

ONE QUARTER OF AN ELEMENTAL FRAME 
Considering member ^ of the frame in Figure 3, the definitive 
differential equation is: 



j. A = a 

and the applicable boundary conditions are: 
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Successive integration of Equation 1, with the application of 



boundary conditions yields: 




fcCT, J, 



Now considering member of the frame in Figure 3, the definitive 
differential equation is: 



Jz - Q C3) 

dXi - 

and the appropriate boundary conditions are: 



j 2 J 3 yJ = 


-H,: 


= -9i 


clx/ IXi-O 




2. 


Jz = 




+ M, 


d XiHxjro 


8 


dYd = 


Qcf 


+ M ,c 




W7, 
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Integration of Equation 3 and the introduction of the above 
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boundary conditions leads to the following expressions: 






\ = Qxt - <9£>z + 
x x 
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From symmetry, the slope at the point X 2 = f/2 is equal to zero. This 
constraint allows the determination of the moment, Mj.. 



M. = 



= q f 2 -p 3 -3cV - C 3 ^)] 

a,# L -P+ cM J 



(5) 



Assuming that the sides of the frame do not differ greatly in length 

or in bending stiffness, the maximum moment occurs at the corners of the 

frame. From Equation 2, the maximum moment can be expressed as: 

1 



K1 



Mtt - &=- 



<3 + 

L 



(6) 



The maximum shear force on the horizontal side of the frame is: 



MAX 



= -Qf 



and the maximum shear force on the vertical side is: 



^MAX ^ 



Equation 6 indicates that J^, J^, f/c, and Q are parameters which 
affect the choice of design from a structural point of view. Based on 
the assumption that the pressure and side dimensions will be dictated 
by requirements for flow and heat transfer, f/c and Q may be fixed for 
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a given design for reasons other than structural. Thus the structural 
designer may have little control over the magnitude of the shear forces. 
The magnitude of the moments may be controlled by varying and J 2 . 
Therefore the structural designer, by carefully choosing the appropriate 
values of member thickness, may be able to optimize the design for the 
idealized frame element* 

In order to study the effect of varying member thicknesses, t^ 
and t 2 > a two dimensional model of the heat exchanger module was in- 
vestigated. A diagram of this model is shown in Figure 4. This model 




FIGURE 4 

TWO DIMENSIONAL HEAT EXCHANGER MODEL 

shows that t^ , t a , t^, and t^ are analogous to the t^ and t 2 shown in 

Figure 3. The effect of varying t , t, , t , and t,- within the model 

a o p i 

chosen is not obvious and cannot be defined in terms of expressions 
derived for moment and shear for the simple frfftie. 

In two dimens ionalizing the model, the top and bottom plates, the 
interior plates, and the closure bars are in reality three dimensional 
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plates which have been idealized as two dimensional beams by assuming 
an elemental cross section of the core. The Plate Flexural Rigidity, 

D, is used for each element to include this effect. 

Slope and deflection methods of analysis are adaptable to matrix 
methods which in turn are easily programmed for solution by the high 
speed digital computer. The computer program, as used, is included 
in Appendix I with a set of answers for the chosen model. This program 
was prepared from considerations of the following. 

The stiffness method of matrix solution is a form of the slope 
and deflection method as proposed by M. L. Pei. j4] The basic element 
for this method consists of a beam under the action of end moments m^ 
and m 2 and end forces Pg and P^ as shown in Figure 5. 
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FIGURE 5 

BASIC ELEMENT, SLOPE AND DEFLECTION ANALYSIS 
Let the beam undergo end rotations ©, and © 2 as well as end displace- 
ments S 3 and . The positive directions of these displacements are 
as defined in Figure 6 . 
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FIGURE 6 

BASIC ELEMENT DEFLECTIONS, SLOPE AND DEFLECTION ANALYSIS 
This beam then, has four degrees of freedom represented by the dis- 
placement matrix: 



Correspondingly, loads are represented by the load matrix: 

rw 

H-k 

l p 4. 

These loads and deflections are related by the stiffness influence co- 
efficient K which depends upon the cross section and the modulus of 
elasticity of the beam. These equations may be expressed conveniently 
in matrix notation as follows : 

M - W {©} 

For the single member example given above, this general equation 
becomes : 
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run, - Kii&i + K»3$3 + ^i4 

rfttz'- Kii&i + Kzz&i + Kja&j + Ki4 &4 

P% * h\ ©» + fat©** K 3 a&i t K34S4 
R} 5 K41 6 » f K41&1 * K 43&3 + |^4 4 S+ 



Hence the "K" matrix is expressed as: 

K 11 Kl 2 Ki3 Ki4- 
K:i Ku 

/(31 Ksi Kn tot 
Kn k43 K 44 

The loading jjYYlj for the structure is known, so the rotations and 
deflections of the member ends, or more generally, the joints in the 
structure, may be determined from: 

{&} = WM 

The application of this method to a general two dimensional frame 
consists of computing the elemental "K M matrix for each element of the 
frame and then adding these elemental "K n matrices to form the total "K" 
matrix for the frame. The "K" matrix is then inverted, as indicated by 
the above equation, to determine the frame joint deflections. 

Using the frame joint deflections the resulting element loads are 
found from : 

[<m} s [k] 

where is the elemental matrix and is the deflection of the 

joints in which the element terminates. These elemental loads are the 
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element internal loads. These internal loads describe the degree of 
relaxation from that of a fixed joint. The external loads or element 
reactions may then be determined from: 

iAftxeb 3 M + Mke Action 

These reactions are the desired results from which stresses may be cal- 
culated. 

The process of calculating loads in this way is very laborious, 
particularly when it is necessary to invert a large "K" matrix. At the 
same time the process involved is purely a mechanical one which may be 
computerized. 
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4. 



RESULTS 



The computer was used to determine the effect of varying t , tjj, 

cL 



and t using the following variations: 
P 

Basic Configuration 
t a = .04 in 
tb = .04 in 
tp = .04 in 



Variations 

.04 in -- .25 in 

.04 in -- .25 in 

.02 in -- .06 in 

Constant 

Constant 

Constant 



tf = .01 in 
f = 0.75 in 
c = 0.5 in 

Maximum normal stresses and shear stresses for the variations shown 
were calculated. The shearing stress variations are functions of fin 
spacing and channel height, neither of which was varied since they were 
assumed to be dictated by thermodynamic considerations. The shearing 
stress variations were small when compared to the normal stress vari- 
ations. For this reason the variations in shearing stresses were assumed 
to be negligible for the parameters being investigated. Graphs showing 
the effect of the changes on shearing stresses are presented in Figures 
8 and 9 . 

Maximum normal stress variations for varying closure bar thickness 
are plotted in Figure 10. Maximum stresses are computed from: 

Maximum normal stress - bending stress plus axial stress. 

As the thickness of the closure bar decreases; 

stress in the closure bar is increased rapidly, 
stress in the top plate is decreased rapidly, 
stresses in the interior plates and fins are decreased 
gradually. 
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Since the stresses in the closure bars increase rapidly, this method 
for reducing stresses in the model is considered marginal. 

Maximum normal stress variations for varying interior plate thick- 
ness are shown in Figure 11. This graph indicates that for an increas- 
ing interior plate thickness, the stresses in all members decrease 
rapidly. However, increasing the thickness of the interior plates 
rapidly reduces the heat transfer effectiveness between channels. Since 
heat transfer is the main purpose of the heat exchanger, this method 
of reducing the stresses must be used judiciously. 

Maximum normal stress variations for varying top plate thickness 
reduces the stresses in all members as long as the top plate is not too 
thick, see Figure 12. Beyond a certain top plate thickness the stresses 
in the top plate and the fins begin to increase. However, further thick- 
ness increase in the top plate is not necessarily detrimental as the closure 
bar stresses continue to fall off rapidly. At the point at which the 
top plate and the fin stresses reach a minimum and the closure bar stress 
reaches a local minimum, there is a significant shift in the location 
of maximum stress both in the closure bars and in the top plate. This 
maximum stress location shift results in a definite dip in the maximum 
stress values. Since variation of the thickness of the top plate has 
no basic undesirable effects, it is most advantageous to make initial 
stress reductions in the model by variation of top plate thickness. 

For any given heat exchanger this computerized method of solution 
is useful in the determination of the range of pressure stresses in 
the module. A computer program and detailed instructions for its use 
are presented in Appendix I. 
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5 . 



CONCLUSIONS AND RECOMMENDATIONS 



The stresses due to internal pressure are a significant part of 
the stresses in a rectangular, compact heat exchanger. The value of 
these stresses for a pressure of 10 atmospheres may be on the order 
of sixty to ninety percent of design allowables* 

The value of stresses in the lattice-like frame which composes the 
heat exchanger core can be controlled through variations in the thick- 
ness of the individual elements involved. 

An optimum ratio of the sum of the stiffness of the top plate 
and interior plates to the stiffness of the closure bars is implied. 

This optimum should consider heat transfer as well as structural re- 
quirements . 

Consideration should be given to increasing the thickness of the 
top plate and/or the thickness of the interior plates in order to 
achieve the required stiffness for these members. This approach can 
lead to an optimum stress level and distribution for an internal pressure 
loading. 

The matrix analysis and the computer program used in this analysis 
should be expanded to provide solutions for three-dimensional problems 
with thermal as well as pressure loads. 

The compatibility of the flange and core is an important structural 
consideration. This should be thoroughly investigated in conjunction 
with further analysis of the core. 
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FIGURE 9. MAXIMUM NORMAL STRESS VERSUS CLOSURE BAR THICKNESS 
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FIGURE 10. MAXIMUM NORMAL STRESS VERSUS INTERIOR PLATE THICKNESS 
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FIGURE 11. MAXIMUM NORMAL STRESS VERSUS TOP PLATE THICKNESS 
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APPENDIX I 



COMPUTER PROGRAM "FRAME" DESCRIPTION AND USE 

Program Frame is a self-contained, Fortran 60, digital computer 
program designed to provide joint deflections and/or rotation and member 
end forces and/or moments for plane rigid frame type structures under 
general loading conditions . 

The geometry of the structure is restricted to a plane structure 
with all members parallel to either the X or Y axis in the usual 
cartesian coordinate system. 

A generalized flow diagram showing the basic mechanics of the 
program follows* 
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The following information is for those who wish to actually use 
the program in its present form. 

Properties such as the modules of elasticity (EMOD) , the area 
moment of inertia (ANER) , and the span may vary from one member to 
another; however, no provision is made for a single member with a 
variable moment of inertia (tapered member), except that an appro- 
priate average might be used if the operator so desires. 

Introduction of actual values of modulus of elasticity and moment 
of inertia will provide results, in appropriate units, of actual de- 
flections and/or rotations. Or, if desired, the relative magnitude of 
modulus of elasticity and moment of inertia, based on selected refer- 
ences may be used instead of using the actual values. In this case, 
the results will indicate the rotation or deflection times a factor 
which is equal to (E re f) x (I re f)* In either case the correct values 
of member end forces and moments in appropriate units will result. 

The required input data is relatively straight forward and simple 
in that the first data card describes the overall structure and sub- 
sequent cards describe the individual members taken one at a time. 
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Care should be taken in compiling the data cards, since seemingly 
insignificant errors will produce entirely erroneous results. 

The program provides analysis for only one structure at a time. 

This structure may have only one complete set of loads applied to it. 

The number of loads comprising this one set, however, is essentially 
unlimited. Provision is made for application of concentrated loads 
and/or distributed loads to any member and in any combination. 

The program output format is self-explanatory except that THETA 
may be actual rotations or deflections if actual values of modulus of 
elasticity and moment of inertia were used, or they may be multiples 
of the reference values if such reference values were used. 

Following is the nomenclature of the program, a copy of the actual 
program in Fortran 60, a diagram for the model actually solved indicating 
member numbers and degrees of freedom (degrees of freedom indicated 
by arrows), data tabulation for the model as shown, a solution for 
deflections and rotations for the model, and the member end moments 
and shears for the model and data cards shown. 

Nomenclature - (Program Proper ) 

AK - Stiffness matrix for an individual member. 

A - Matrix used for obtaining the stiffness matrix for the 
complete structure . 

C - Stiffness matrix for the complete structure. 

THETA - Matrix of deflections and/or rotations. 

AKS - Matrix used to store individual AK matrices and M, N, 

K, L, and end moments and forces due to external loads. 

DELTA - THETA matrix for a single member. 

F - Matrix of end forces and moments for a single member. 
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Required Input Data 

Data Card #1 Contains values of NX, NM, EP, punched in the 
format (2110, E10.4). 

Data Card #2 Describes the first member to be considered. 

It contains values of M, N, K, L, EMOD, ANER, 
SPAN, punched in the format (412, 2E10.5, 

1F10.5) . 

Data Card #3 Defines the number of loads (NL) applied to 

the member described in card #2. NL is punched 
in the format (12). If no load is applied to 
the member being considered, a zero must be 
punched in column number 2. 

Data Card #4 Contains values of KFM, W, SPA 1, SPA 2, punched 
in the format (12, 3F10.5). If data card #3 
set NL equal to zero, this card (#4) must be 
omitted. The total number of #4 cards used at 
this point must equal NL from card #3. 

Subsequent Data Cards The input data card sequence is now 
repeated starting with card # 2 above. Each 
member, in turn, must be completely described 
with its properties and loads before proceeding 
to the next member. 

Nomenclature - (Input Data ) 

NX - Number of degrees of kinematic deficiency. 

NM - Number of members or spans comprising the structure. 

EP - Error parameter, minimum acceptable value of the deter- 
minant of the stiffness matrix. If determinant is less 
than EP the program will stop and print M matrix singular." 



35 



M - The number assigned to the rotational freedom at the left 
end of a member. 

N - The number assigned to the rotational freedom at the right 
end of a member. 

K - The number assigned to the translational freedom at the 
left end of a member. 

L - The number assigned to the translational freedom at the 
right end of a member. 

EMOD - Modulus of elasticity for a member. May also be used as 
a coefficient of some reference modulus. 

ANER - Area moment of inertia for a member. May also be used as 
a coefficient of some reference moment of inertia. 

SPAN - The length or span for a member. 

NL - The number of external loads (concentrated and/or uniform) 
applied to the single member being considered. 

KFM - Kind of load applied (1 = concentrated; 2 = uniform). 

W - Magnitude of the external load. 

SPA 1 - Distance from left end of member to the location of a 
concentrated load . 

SPA 2 - Distance from a concentrated load to the right end of a 
member . 

Basic Element 




Note: Vertical members are rotated 90 degrees clockwise to determine the 

orientation indicated above. 
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The structure is sketched and all degrees of kinematic deficiency 
are numbered. 

The individual members are viewed one at a time, as shown on 
previous page, and M, N, K, and L for each member are set equal to the 
number which was used to designate the corresponding degree of freedom 
for that member. 

Load Parameters 



KFM = 1 
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PROGRAM FRAME 
C RUN 110 

READ lit NX *NM* EP 
11 FORMAT (2110, E10.A) 

DIMENSION AK ( A » A ) , A ( 56 » 56 ) ,C ( 56 ♦ 56 ) » THETA ( 56 ) , AKS ( 56 ♦ 2A ) 
1 DELTA ( A ) ♦ F ( A ) 

DO 90 10 * 1» NM 

READ 20* Mt N » K* L, EMOD. ANER » SPAN 
20 FORMAT (AI2. 2E10*5, 1F10.5) 

DO 23 1*1,56 
DO 23 J* 1*56 
23 A ( I * J ) = 0. 

DO 26 1*1, A 
DO 26 J*l, A 
26 AK ( I * J ) * 0. 

AK ( 1 » 1 ) * ( A • *EMOD*ANER ) /SPAN 
AK ( 2 *2 ) *AK ( 1*1) 

AK ( 1 » 2 ) = ( 2.*EM0D*ANER)/SPAN 
AK ( 1 *3 ) = (6,*EMOD*ANER)/(SPAN**2) 

AK ( 2 * 3 ) * AK ( 1 ♦ 3 ) 

AK ( 1 ♦ A ) = -AK (1*3) 

AK ( 2 ♦ A ) * -AK ( 1,3) 

AK ( 3 * 3 ) = ( 12,*EMOD*ANER)/(SPAN**3) 

AK( A,A ) = AK ( 3 , 3 ) 

AK ( 3 » A ) = -AK (3*3) 

AK (2*1) =AK ( 1*2) 

AK(3»1)*AK( 1*3) 

AK(3*2)=AK(2»3) 

AK ( A , 1 ) =AK ( 1 » A ) 

AK (A,2)=AK( 2»A) 

AK ( A , 3 ) = AK ( 3 , A ) 
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AO 

A2 

A3 

A5 

A6 

AS 

A9 



AKS ( 10,17) =M 
AKS ( I 0 , 1 8 ) =N 
AKS ( 10* 19 ) =K 
AKS ( 10* 20 ) *L 
KK* 1 

DO 3615 LL= 1 * A 
DO 3615 J J= 1 » A 
AKS( IO,KK)=AK(LL»JJ) 
KK=KK+1 

IF (M) 103»A8»38 
A ( M *M ) * AK (1*1) 

IF (N) v 103* A2» AO 
A ( M ,N ) = AK ( 1 *2 ) 

A ( N ,M ) * AK (1*2) 

IF (K) 103* A5, A3 
A ( M*K) = AK (1*3) 

A ( K *M) * AK ( 1 ♦ 3 ) 

IF (L) 103, A8 * A6 
A ( M*L ) * AK ( 1 , A ) 
A(L*M) * AK ( 1 » A ) 

IF (N) 103, 56* A9 
A(N*N) * AK (2*2) 
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IF (K) 103 * 53 ♦ 51 
51 A ( N ♦ K ) = AK ( 2 » 3 ) 

A ( < * N ) = AK ( 2 » 3 ) 

53 IF (L) 103 » 56, 54 

54 A ( N , L ) * AK ( 2 ,4 ) 

A ( L ,N ) * AK ( 2 ,4 ) 

56 IF (K) 103, 61, 57 

57 A ( K , K ) = AK ( 3 ,3 ) 

IF (L) 103, 61, 59 

59 A ( K ,L ) = AK (3,4) 

A ( L ♦ K ) = AK ( 3 ,4 ) 

61 IF < L ) 103, 63, 62 

62 A ( L ,L ) = AK ( 4 ,4 ) 

63 DO 65 1=1,56 
DO 65 J=1 ,56 

65 C( I , J) = C( I , J) + A ( I , J) 

NPM=NX +1 
READ 681, NL 
681 FORMAT (121 

IF (NL) 103, 90, 69 
69 DO 90 KO = 1,NL 

READ 71, KFM ♦ W, SPA1, SPA2 
71 FORMAT (112, 3F10.5) 

IF (KFM-2) 73, 82, 103 

73 IF ( M ) 103 ♦ 741 , 74 

74 C ( M , NPM ) = (W*SPA1*SPA2**2)/(SPAN**2) + C(M,NPM) 

741 AKS (I0»21)=-(W*SPA1*SPA2**2)/(SPAN**2) +AKS( 10,21) 

IF (N) 103,761 ,76 

76 C ( N , NPM ) = -(W*SPA2*SPAl**2)/( SPAN**2 ) + C(N,NPM) 

761 AKS(I0,22)=(W*SPA2*SPA1**2)/(SPAN**2) +AKS(IO,22) 

IF (K) 103,781,78 

78 C ( K , NPM ) = ( (W*SPA2**2 )*( 3.*SPAl + SPA2 ) ) / ( SPAN** 3 ) + C(K»NPM) 
781 AKS ( 10,23) =- ( ( W*SPA2**2 ) * ( 3 . *SPA 1+SPA2 ) )/(SPAN**3) +AKS( 10,23) 
IF (L) 103,801,80 

80 C ( L ♦ NPM ) = ( (W*SPA1**2 )*( 3.*SPA2 + SPA1) )/(SPAN**3) + C(L,NPM) 
801 AKS (10, 24 ) = -*( (W*SPA1**2)*(3.*SPA2+SPA1) )/(SPAN**3) +AKS( 10,24) 
GO TO 90 

82 IF (M) 103,831,83 

83 C ( M ♦ NPM ) = ( W* SPAN** 2 ) / ( 12 * ) + C(M,NPM) 

831 AKS(IO,21)=-(W*SPAN**2)/(12.) +AKS(IO,21) 

IF (N) 103,851,85 

85 C ( N , NPM ) = -(W*SPAN**2 )/( 12. ) + C(N,NPM) 

851 AKS( 10,22 )= ( W*SPAN*#2 ) / ( 12 • ) +AKS(IO,22) 

IF(K)103,871,87 

87 C ( K , NPM ) = (W*SPAN)/(2. ) + C(K»NPM) 

871 AKS( IO,23)=-( W*SPAN ) / ( 2 • ) +AKS(IO»23) 

IF ( L ) 103,891,39 

89 C ( L , NPM ) = ( W#SPAN ) / ( 2 • ) ■ + C(l»NPM) 

891 AKS( IO»24)=-( W*SPAN ) / ( 2 • ) +AKS(IO»24) 

90 CONTINUE 
NP * 1 

CALL GAU992 !NX» NP» EP» €» THETA, Kl) 

GO TO (93, 99) , Kl 
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93 

94 

321 

96 

97 

9072 

9074 

730 

731 

976 

979 

9711 

9712 

9714 

9716 

9717 

9719 

9721 

9722 

9724 

9726 

9727 

9730 

9735 

9736 

99 



PRINT 94 

FORMAT ( ////44H1N0TE 1. THETA IS JOINT ROTATION/DEFLECTION/) 

PRINT 321 

FORMAT ( 29H NOTE 2. ANSWER = El X THETA////) 

DO 96 K=1 »NX 

PRINT 97, K ♦ THETA ( K ) 

FORMAT (11H0EI X THETA 13, 3H = F20.15) 

PRINT 9072 

FORMAT ( ////43H1MEMBER END MOMENTS AND SHEARS APPEAR BELOW/) 
PRINT 9074 

FORMAT ( 47H CLOCKWISE MOMENTS AND DOWN SHEARS ARE POSITIVE////) 
PRINT 730 

FORMAT! ////43H LEFT RIGHT LEFT RIGHT) 

PRINT 731 

FORMAT (43H MEMBER MOMENT MOMENT SHEAR SHEAR//) 

DO 9735 I M= 1 * NM 
KJ=> 1 

DO 976 LL= 1 *4 
DO 976 J J = 1 » 4 
AK ( LL ♦ J J ) = AKS ( IM.KJ) 

KJ = K J +1 



M=AKS ( I M » 1 7 ) 

IF (M) 103.9711.979 
DELTA! 1 ) =THETA ( M ) 

GO TO 9712 
DELTA! 1 )=0. 

N= AKS ( I M ♦ 1 8 ) 

I F ( N ) 103,9716.9714 
DELTA! 2 ) = THETA ( N ) 

GO TO 9717 
DELTA! 2 )=0. 

K = AKS ( I M » 1 9 ) 

IF (K) 103,9721,9719 

DELTA ( 3 ) =THETA ( K ) 

GO TO 9722 
DELTA! 3 )=0. 

L=AKS ( I M ♦ 20 ) 

IF (L) 103,9726,9724 
DELTA(4)=THETA(L) 

GO TO 9727 
DELTA ( 4 ) =0 • 

DO 9730 KB= 1,4 
F(KB)«0. 



DO 9730 JB= 1,4 
F ( KB ) =F ( KB ) +AK ( KB , JB ) * 
F ( 1 ) =F ( 1 ) +AKS ( I M , 2 1 ) 



F ( 2 ) =F ( 2 ) + 
F ( 3 ) =F ( 3 ) + 
F ( 4 ) =F ( 4 ) + 
PRINT 9736. 
FORMAT (14, 
GO TO 105 
PRINT 100 



AKS! 


IM, 


22) 


AKS 


( IM 


,23 ) 


AKS! 


IM. 


24) 


I M ♦ F 


( 1 ) 


♦ F ( 2 


F9.3 


♦3F7.3) 



DELTA! JB) 



♦ F ( 3 ) , F ( 4 ) 
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100 FORMAT (16H MATRIX SINGULAR) 

GO TO 105 

103 PRINT 104 

104 FORMAT (12H INPUT ERROR) 

105 CONTINUE 
STOP 
END 

SUBROUTINE GAUSS2 ( N »M» EP » A , X » KER ) 

DIMENSION A ( 56 ♦ 56 ) * X ( 56) 

NPM=N+M 

10 DO 34 L=1»N 
KP = 0 

Z = 0.0 

DO 12 K=L ♦ N 

IF ( Z-ABSF ( A( K»L) )) 11 *12* 12 

11 Z=ABSF(A(K»L) ) 

KP = K 

12 CONTINUE 

I F ( L-KP) 13*20t20 

13 DO 14 J = L * NPM 
Z=A(L» J) 

A(L,J)=A(KP*J) 

14 A ( KP » J ) =Z 

20 I F < ABSF(A(L»L) ) -EP ) 50 » 50 » 30 

30 IF(L-N)31» 40*40 

31 LP1=L+1 

DO 34 K=LP1*N 
I F ( AC K *L ) )32.34,32 

32 RATIO=A(K»L)/A(LtL) 

DO 33 J=LP 1 »NPM 

33 A(K*J)=ACK»J) -RAT 1 0 # A ( L ♦ J ) 

34 CONTINUE 

40 DO 43 I=1»N 
I I =N+ 1- 1 

DO 43 J= 1 » M 
JPN=J+N 
S = 0.0 

I F ( I I-N)41*43»43 

41 IIP1-II+1 

DO 42 K= 1 1 PI *N 

42 S = S+A ( I I * K )*X ( K) 

43 X < 1 1 ) = C A ( I I * JPN ) -S ) / A ( I I * I I ) 

KER = 1 

RETURN 
50 KER = 2 
END 
END 
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HEAT EXCHANGER MODEL FOR COMPUTER SOLUTION 
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NOTE 1 — T H ETA - TS~ J 01 NT - ROT A T I ON / DE FLEC'T ION 

NOTE 2. ANSWER = El X THETA 



El X THETA 1 = -.059757168141914 

-El— X— TH ETA 2 = -“. 039902351722048“ 

E I X THETA 3 = -.000000000003752 

El X THETA 4 = .039902351729324 

El X THETA 5 = .059757168116448 

"E'l X TH2i5A 6 = .020913817393648- 

E I X THETA 7 = -.084237065282650 

El X THETA 8 = -.000181155194702 

El X THETA 9 = -.030079360703894 

El— X-THETA— 10-= T00000000000 682 1“ 

EI X THETA 11 = .030079360741638 

El X THETA 12 = .000181155228802 

El X THETA 13 = .084237065302659 

-£7 X— TH ETA— 1-4— = -“.“02 0913817406391- 

EI X THETA 15 = -.023057530810092 

El X THETA 16 = .038969665726654 

El X THETA 17 = -.041627451475506 

ET—X— THETA— 1 8~= .0 1238007425 331 S~ 

E I X THETA 19 = -.000000000023057 

El X THETA 20 = -.012380074174139 

El X THETA 21 = .041627451528257' 

-El— X— TH E T A— 22 = =■ ."0 38969 665821242“ 

EI X THETA 23 = .023057530988808 

El X THETA 24 = .025138163619431 

El X THETA 25 = .003051 132390510 

-£ I— X— TH E T A - 26 = .05 3 1 46 5381 0723 5“ 

EI X THETA 27 = -.011316232965555 

El X THETA 28 = - .000000000938599 

El X THETA 29 = .011316233752268 

-£I X THET A~ 30~ = -.053146539722547“ 

E I X THETA 31 = -.003051131717825 

El X THETA 32 = -.025138164648979 

El X THETA 33 = .049425681714638 

-£-1— X— THETA-3 4— = . 04 77 1 302039 534 9“ 

E I X THETA 35 = .018272389301728 

El X THETA 36 = -.018272389324466 

El X THETA 37 = -.047713020370793 

-E I— X— TH ETA 38 -.04 942568167916 8“ 

E I X THETA 39 = -.039716653672258 

El X THETA 40 = -.054848966682594 

El X THETA 41 = -.039716653668620 

-Ei— X— TH E T-A-4-2-= -“.00 5222375 323 56 9“ 

E I X THETA 43 = -.019661487614485 

EI X THETA 44 = -.019661487617668 

EI X THETA 45 = -.005222375324706 

EI— X— T H E T A - 46 = . 011389984892958“ 

E I X THETA 47 = .012368798751140 

EI X THETA 48 = .011389984899324 

EI X THETA 49 = .018564758140656 

-Ei— X— TH ETA— 5^0 — = rO 4-4-1 99 61 4962963“ 

E I X THETA 51 = .044199614952049 

EI X THETA 52 = .018564758128377 
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MEMBER EN B— MOMENT-S— AND-SH EA RS~A P P E A R— BELOW 

CLOCKWISE MOMENTS AND DOWN SHEARS ARE POSITIVE 



MEMBER 



LEFT 

MOMENT 



RIGHT 

MOMENT 



LEFT 

SHEAR 



RIGHT 

SHEAR 



1 

2 

' 3 
- 4 - 



5 

6 
7 

- 8 - 



9 

10 
1 1 
- 1 - 2 - 
1 3 

14 

15 
9 6 - 

17 

18 
1 9 

- 20 - 



21 

22 

25 

- 2 - 4 - 



25 

26 
27 
- 2 - 3 - 



29 

30 

31 
- 3 - 2 - 

33 

34 

35 
- 36 - 



37 

38 

39 
- 4 - e - 

41 

42 

43 
- 44 - 

45 

46 

47 
- 4 - 6 - 



49 

50 

51 



6.28 

- 19.29 

- 21.13 

-+ 9 .- 33 - 



-.03 

4.58 
- 8 . 14 

9 .- 44 - 

- 10.73 
9.42 
- 8 . 16 

4 -.- 5 - 1 - 

1.48 

- 4.12 

7.57 

9v88 - 

10.40 

- 9.88 

7.56 

4 -.- 1 - 2 - 

- 4.06 
4.98 
- 7.15 
- 9 -.- 90 - 



- 10.38 

9.90 
- 7 . 11 

5 .- 00 - 

- 7.09 

8.20 

17.92 

-+ 7 -.- 93 - 



8.25 
2.55 
- . 02 

. 00 - 

.02 
- 2.55 
2 . 32 

.- 00 - 

-.00 

.00 

.00 

- 2S - 32 - 



4.00 

.03 

-.00 

— -.- 03 - 
- 4.00 
7.09 
.05 



19.33 
21 . 13 
19.29 
— 6 .- 28 - 
- 4.51 
8.16 
- 9.42 
- 10 .- 73 - 
- 9.44 
8.14 
- 4.58 
. 03 - 



4.12 
- 7.56 
9.88 
HO -.- 40 - 
9.88 
- 7.57 
4.12 
— 1 -. 48 - 
- 5.00 
7.11 
- 9.90 
- 10 -.- 38 - 
- 9.90 
7.15 
- 4.98 
—4 . 06 - 
- 8.25 
- 17.93 
- 17.92 
— 8 . 20 - 
7.09 
- 6.28 
-.04 
. 00 - 



.04 

6.28 

- 2.52 

. 06 - 

-.02 

.02 

.06 

- 2 -.- 52 - 



- 3.80 

-.02 

-.00 

.-02- 

3.80 
.07 
.03 



90.39 

58.71 

53.79 
— 22 .- 11 - 

- 12.1 1 
33.98 
- 46.84 
— S3 -. 79 - 
- 53.79 
46.84 
- 33.98 
— T2 -. il - 
14.94 
- 31.15 
46.54 
— 54 . 08 - 
54.08 
- 46.54 
31.15 
— 14 . 94 - 
- 24.1 7 
32.24 
- 45.45 
— 54 . 08 - 
- 54.08 
45.45 
- 32.24 
— 24 . 17 - 
- 69.05 
- 69.22 
- 56.25 
— 43 .; 

12.80 
30.04 

-.12 
. 00 - 



.12 

- 30.04 

37.09 

. 14 - 

-.05 
.05 
.14 
— 37 . 09 - 
37.90 
.03 
-.00 
-. 03 - 



- 37.90 
51 . 82 - 
.15 
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- 54 . 08 - 
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APPENDIX II 



OTHER METHODS OF ANALYSIS 



Prior to the decision to use the two-dimensional frame analysis 
for the internal pressure stresses a number of other methods were in- 
vestigated. Generally, the approach was one of searching for the correct 
design philosophy and attempting to adapt simple solutions to the prob- 
lem at hand. 

The methods outlined herein seemed to apply to the problem but in 
each case the results became too complex to provide useful design in- 
formation or the results required specific knowledge of some particular 
parameter which was unavailable. 



The first approach was to consider the top plate of the heat ex- 
changer to be fixed on the two sides adjoining the side closure bars, 
to be resting on an elastic foundation, with modulus K lb/in, along its 
span, and to be loaded with a uniform pressure load, Q lb/in. The elastic 
foundation was assumed to be the result of the arrangement of interior 
plates and fins. Assuming no change in geometry or other parameters 
with length, the plate was analyzed as a beam considering the stiffness 
to be El/(1-V^) instead of simply El. 

The expression for the moment at the end of this beam is: 



This result shows that the modulus of the elastic foundation, K, is of 
primary importance. In truth, this modulus is a function of such para- 
meters as fin spacing, width of the heat exchanger, and thickness or 



TOP PLATE ON AN ELASTIC FOUNDATION 




where 
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stiffness of the interior plates. No adequate analytical means of 
estimating the magnitude of the foundation modulus was found. It is 
recommended that K be determined by testing actual modules. 

DISCONTINUITIES BETWEEN HEADER AND CORE 
A number of attempts were made to determine the value of discon- 
tinuity stresses at the joint between the header or flange and core. 

In general, this method attempted to extend the classical solutions for 
cylindrical pressure vessels to include rectangular pressure vessels. 

This approach might be possible if the corners of the rectangular heat 
exchanger could be rounded off slightly such that a continuous mathematical 
function could be used to define the surface of the heat exchanger. The 
best approach with the square corners appears to be with the use of an 
elemental lengthwise strip taken at the location of the greatest dis- 
continuity in deformation between the core and the flange as shown below. 




If Q q and Mq could be determined, the elemental strip of the top 
plate of the core could be analyzed as a beam on an elastic foundation 
with end loads Q q and M Q . Difficulties, in connection with finding Q q 
and M 0 , the magnitude of the elastic foundation modulus, whether the 
elemental strip of the core could be considered semi-infinite in length, 
and the correct method for determining the between the header and 
core, were encountered. 
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FRAME ANALYSIS 



A number of lattice -type frames were analyzed and the analysis for 
the simple frame with two axis of symmetry did, in fact, provide the 
impetus for continuing with the complex frame analysis as used. This 
analysis provided fundamental information as to the effect of the various 
parameters involved in the core. 

Frames with various numbers of non-yielding tie bars, which repre- 
sented the interior plates, were analyzed. The solutions neglected the 
effect of fins. This approach was abandoned because of the inability 
to adequately account for the effect of the fins. 

SANDWICH CONSTRUCTION 

Considerable effort was expended in researching publications dealing 
with sandwich construction. Much is available on the subject in the 
literature but no reference was found which dealt with sandwich or honey- 
comb structures loaded with an internal pressure. 
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